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Abstract 

In deeply inelastic lepton-nucleus scattering (DIS), the average jet transverse 
momentum is broadened because of multiple scattering in the nuclear medium. 
The size of jet broadening is proportional to the multi-parton correlation 
functions inside nuclei. We show that at the leading order, jet broadening 
in DIS and nulcear enhancement in di-jet momentum imbalance and Drell- 
Yan (qj,) share the same four-parton correlation functions. We argue that jet 
broadening in DIS provide an indenpendent measurement of the four-parton 
correlation functions and a test of QCD treatment of multiple scattering. 

11.80.La, 12.38.Bx, 13.85.Qk, 24.80.-x 



Typeset using REVTgX 



*This work was partially supported by the U.S. Department of Energy under contract No. DE- 
FG02-93ER40764. 

1 



I. INTRODUCTION 



When a parton propagates through a nuclear matter, the average transverse momentum 
may be broadened because of multiple scattering of the parton. Nuclear dependence of 
such broadening provides an excellent probe to study QCD dynamics beyond the single 
hard-scattering picture. Reliable calculation of multiple scattering in QCD perturbation 
theory requires to extend the factorization theorem JO beyond the leading power. Qiu and 
Sterman showed that the factorization theorem for hadron-hadron scattering holds at the 
first-nonleading power in momentum transfer O , which is enough to study double scattering 
processes in hadronic collisions. According to this generalized factorization theorem, the 
double scattering contribution can be expressed in terms of universal four-parton (twist- 
four) correlation functions. The predictions of the multiple scattering effects rely on the 
accurate information of the four-parton correlation functions. Various estimate of the size 
of these four-parton correlation functions were derived recently [0-0 • It is the purpose of this 
paper to show that jet broadening in deeply inelastic lepton-nucleus scattering (DIS) can be 
used as an independent process to test the previous estimate of the four-parton correlation 
functions. 

The four-parton (twist-4) correlation functions are as fundamental as the normal twist-2 
parton distributions. While twist-2 parton distributions have the interpretation of the prob- 
ability distributions to find a parton within a hadron, four-parton correlation functions pro- 
vide information on quantum correlations of multi-partons inside a hadron. Just like normal 
twist-2 parton distributions, multi-parton correlation functions are non-perturbative. QCD 
perturbation theory cannot provide the absolute prediction of these correlation functions. 
But, due to factorization theorems, these correlation functions are universal; and can be 
measured in some processes and be tested in other processes. 

Luo, Qiu and Sterman (LQS) calculated the nuclear dependence of di-jet momentum 
imbalance in photon-nucleus collision in terms of nuclear four-parton correlation functions 
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X \ <Pa\MO) 7+ F/fe ) ^ + °(»r) lMiT)M ; (la) 



and 



x 4r <P-4|F„ + (0) F„+(jfe) F+'fef) F+"(,r)M ■ (lb) 

By comparing the operator definitions of these four-parton correlation functions and the 
definitions of the normal twist-2 parton distributions, LQS proposed the following model 

T f/A (x) = X 2 A^<f> f/N (x) , (2) 

where 0//tv(x) with / = q,q,g are the normal twist-2 parton distribution of a nucleon, and 
A is a free parameter to be fixed by experimental data. Using the Fermilab E683 data on 
di-jet momentum imbalance, LQS estimated the size of the relevant four-parton correlation 
functions to be of the order A 2 « 0.05 ~ 0.1 GeV 2 §. 

The nuclear enhancement of the average Drell-Yan transverse momentum, A(gfA, also 
depends on the similar four-parton correlation functions || . However, Fermilab and CERN 
data on nuclear dependence of the Drell-Yan A(gfi) prefer a much smaller size of the four- 
parton correlation functions 0,[7]||. Actually, as we will show below, the Drell-Yan data 
favor the four-parton correlation functions about five times smaller than what was extracted 
from the di-jet data. This discrepancy may result from different higher order contribution 
to the Drell-Yan A(gfi) and the di-jet momentum imbalance. We have pure initial-state 
multiple scattering for the Drell-Yan process, while pure final state multiple scattering for 
a di-jet system. It is necessary to study the high order corrections to these two observables 
in order to test QCD treatment of multiple scattering. It was argued in Ref. P,|10f that 



the high order contribution to the Drell-Yan nuclear enhancement is important. Meanwhile, 
it is also important and necessary to find different observables that depend on the same 
four-parton correlation functions. 



The jet broadening in deeply inelastic lepton-nucleus scattering provides an independent 
test of the size of the four-parton correlation functions. At the leading order in deeply 
inelastic scattering, a quark scattered of the virtual photon forms a jet (known as the 
current jet) in the final state. Because of the multiple interaction when the scattered parton 
propagates through the nuclear matter, the transverse momentum of the final state jet is 
broadened in deeply inelastic lepton-nucleus scattering |l]J . We show below that the size of 
the jet broadening is directly proportional to the four-parton correlation functions. These 
four-parton correlation functions are the same as those appeared in the nuclear enhancement 
of the di-jet momentum imbalance. Current data for the di-jet momentum imbalance and 
the Drell-Yan A(gf,) provide two different sizes of the four-parton correlation functions. 
Using these data, without additional parameters, we predict a range of the jet broadening 
in DIS. Existing data from Fermilab E665 and future HERA nuclear beam data can provide 
a direct test of QCD treatment of multiple scattering. 

The rest of this paper is organized as follows. In Sec. II, we derive the nuclear enhance- 
ment of the Drell-Yan (gf,) by using the same method used by LQS to derive the di-jet 
momentum imbalance, and show that our result is consistent with that presented in Ref . [|J . 
In Sec. Ill, using the same method, we derive jet broadening in DIS in terms of the same 
four-parton correlation functions. Finally, in Sec. IV, we use the values of A 2 extracted from 
data on di-jet momentum imbalance and Drell-Yan A(q^) to estimate jet broadening in DIS. 
We also present our discussions and conclusions in this section. 



II. NUCLEAR ENHANCEMENT FOR DRELL-YAN (Q 2 T ) 

Consider the Drell-Yan process in hadron-nucleus collisions, h(p') +A(p) — > £ + £~(q) +X, 
where q is the four-momentum for the virtual photon 7* which decays into the lepton pair. 
p' is the momentum for the incoming beam hadron and p is the momentum per nucleon for 
the nucleus with the atomic number A. 

Let qx be the transverse momentum of the Drell-Yan pair, we define the averaged trans- 
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verse momentum square as 

In Eq. (^|), Q is the total invariant mass of the lepton pair with Q 2 = q 2 . The transverse 
momentum spectrum, da/dQ 2 dq^, is sensitive to the multiple scattering, and has the A 1 ^ 3 - 
type nuclear size effect. If we write the cross section as 

<?hA = °hA + °hA + ••• > ( 4 ) 

where a^ A , <j® a , and represent the single, double and higher multiple scattering, re- 
spectively. The single scattering is localized, hence, a^ A does not have a large dependence 
on the nuclear size. Therefore, a^ A ~ AahN with N represents a nucleon. As we will demon- 
strate explicitly, the inclusive cross section da/dQ 2 has much weaker nuclear dependence 
after integration over the whole momentum spectrum. Hence, 



J dq 2 T ■ q 2 T 



^ a hA / d^hA ^ / n 2\hN /r\ 

dQ 2 dq 2 I dQ 2 ~ {qT) • [b) 
To extract effect due to multiple scattering we introduce the nuclear enhancement of the 
Drell-Yan (gf,) as 

A(g|) = (q 2 T ) hA - (q 2 T ) hN • (6) 

If we keep only double scattering and neglect contribution from the higher multiple scatter- 
ing, we have 

From our definition, A(gf,) represents a measurement of QCD dynamics beyond the tradi- 
tional single-hard scattering picture. The nuclear enhancement of Drell-Yan (qj), defined in 
Eq. (|6|), is a result of multiple scattering between the incoming beam parton and the nuclear 
matter before the Drell-Yan pair is produced. 

Since we are interested in the averaged transverse momentum square, we do not need 
to know the angular information of the Drell-Yan lepton pair. After integration over the 
lepton's angular information, the Drell-Yan cross section can be written as 



da, 



2a P 



hA^t+t- 



3g 2 



(8) 



where is the hadronic tensor []E|. To simplify the notation of following derivation, we 
introduce 



W hA ^(q) = -g^W^Jq) 



At the leading order, the single scattering contribution to WhA->-y* (<?) is given by 



(9) 



W^ r (q) =J2 J dx' (j)q /h {x') J dx(j) q/A {x) 



2 27TQf e 77T, 



(10) 



where superscript "S" signals the single scattering; and the corresponding inclusive single 



scattering cross section is 12 



da 



~dq r ~ 



a oJ2 e lJ dx' (f)q/ h (x') J dx(j) q/A {x) 5(Q 2 - xx's) 



(11) 



with s = (p + p') 2 and the Born cross section 



Aira 2 

em 



9g 2 



(12) 



Consider only double scattering inside the nuclear target, we can factorize the double 
scattering contribution to W(q) as 



Ka^(i) = E / dx'<j) f/h {x') ■ Wf A ^{x\q) 



(13) 



where superscript "D" indicates the double scattering contribution, and / sums over all par- 
ton flavors. In Eq. (|13|), Wf A ^* (x', D) is the hadronic contribution from double scattering 
between a parton / and the nucleus. At the leading order in a s , Wj^^* (x 1 , D) is given 
by the Feynman diagrams in Fig. [|. According to the generalized factorization theorem , 
Wf A _ yy *(x' , D) can be factorized as 



(14) 
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where T D is the phase space factor for the double scattering. For diagram shown in Fig. [l]a, 
for example, the phase space factor is given by 

= • (2?r) + x lP + x'p' + k T -q) 



(27T) 

■\- 5(x + Xl - ^ - 5{q 2 T - k 2 T ) dQ 2 dq 2 T . (15) 



X'S X'S X'S 



At this stage of the derivation, T and H in Eq. flU]) both depend on gauge choice, even 
though W{q) D is gauge invariant. In this paper, Feynman gauge is used in our calculation. 
In Eq. (|H|), the hadronic matrix element 

rpmt , x Ml/" ^1" dy£ d 2 y T 

T [ >{x,x u x 2 ,k T ) 



2tt 2tt 2tt (2tt) 2 

x g«P + S/" e ixip + fa -J/2 ) gia^p+J/j" e ik T -y T 

x^{p A \A + {y^,0 T )^ q {0) 1 + ^ q {y-)A + {yY,y 1T )\PA) • (16) 

Because of the exponential factors in Eq. (0), the position space integration, dy~'s cannot 
give large dependence on nuclear size unless the parton momentum fraction in one of the 
exponentials vanishes. If the exponential vanishes, the corresponding position space integra- 
tion can be extended to the size of whole nucleus. Therefore, in order to get large nuclear 
enhancement or jet broadening, we need to consider only Feynman diagrams that can pro- 
vide poles which set parton momentum fractions on the exponentials to be zero ||. At 
the leading order, only diagrams shown in Fig. [I] have the necessary poles. These diagrams 



contribute to the double scattering partonic part H{x, xi,x 2 , kr, x'p',p, q) in Eq. (|T4|). 

For the leading order diagrams shown in Fig. [I], the corresponding partonic parts have 
two possible poles. For example, for diagram shown in Fig. [I], the partonic part has following 
general structure 

{x'p' + X\p + k T ) 2 + it {x'p' + {x\ — x 2 )p + k T ) 2 — ie ^ 
« P l — » ( 17 ) 

where H^ u is the numerator which is proportional to a quark spinor trace. 



The first step to evaluate the Eq. ([14]) is to carry out the integration over parton momen- 
tum fractions dxdx\ dx 2 . Using one of the 5-functions in phase space factor in Eq. flT5| ) to 
fix dx integration, and the two poles in the partonic part in Eq. (|T7| ) to perform the contour 
integration for dx\ dx2, we obtain 

I a = Jdxd Xl dx 2 e ^p + y- e^ivi-v^) e ^p + y^ 

1 1 . h\ Q 2 . 

X U • Uz 0{X + Xi ; — J 

X\ r + ie X\ — X2 j — ^e x 6 x b 

1 x's 1 * x's 

= (4tt 2 ) 6(y- - yr) 9(-y 2 ) e l(T/x ' )p+y ~ e H^/^ s )) P +(y^- y -) ? ^ 

with t = Q 2 /s. For the diagram shown in Fig. [l]b, we have a slightly different phase space 
factor 

rf = \- S(x + %2 - S(q 2 T ) dQ 2 dq 2 T . (19) 
x's x's 

The corresponding partonic part H b has slightly different poles 

H b oc 1 l — . (20) 

1 x's 

Integrating over parton momentum fractions, dxdx\ dx 2 , we have 

h = J dx dx x dx 2 e lxp+y ~ e ixip+iy ^- y ^ e lX2P+y * 
1 1 , ( ^ Q\ 

X 72 • 0{X + X 2 — J 

x ^ + ie x 2 + ie x's 

1 x's 

= (-4tt 2 ) 9(y 2 - yl) 6{y~ - y 2 ) e l(r/x ' )p+y ~ e i^ T l{x's)) P +{y-- y -) _ ^ 
Similarly, for diagram shown in Fig [Tjc, we have 

I c = (-An 2 ) %f - y 2 ) 9(-y^) e l{r/x ' )p+y ~ e ^ 2 T /{x's)) P +(y-- y ^) _ ^ 



It is easy to show that the numerator spinor trace gives the same factor for all three 
diagrams. Therefore, total contribution to jet broadening is proportional to I a + h + h, and 
this sum has the following feature 
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d(T hA 

dQ 2 dq\ 



oc I a + I b + I c 



oc 0{y- 



+ 



o(y- 



Vx ) 
-Vx 



-Ui ) 

~V2 



-9(y--y-)9(y- 



V2 



(23) 



Ofa-VM-Hl) *(<&) ■ (24) 



It is clear from Eq. (|2~3| ) that for the inclusive Drell-Yan cross section, da/dQ 2 , the double 
scattering contribution, da® A /dQ 2 , does not have large dependence on the nuclear size. The 
integration over dq\ eliminates the first term in Eq. fl24f ), while the second term is localized in 
space if t/x' is not too small. When the t/x' is finite, and p + is large, the exp[i (r/x')p + y~] 
effectively restricts y~ ~ 1/((t/x')p + ) — > 0. When y~ — > 0, the combination of three pairs 
of ^-functions in Eq. fl24|) vanishes, 



oc 







7 1 '-2 dg 2 

%~ - vx) H-V2) - % 2 ~ - vx) %~ - V2) - %r - yj) 0(-2/D 

as y — ► . 



(25) 

Physically, Eq. fl2"5D says that all integrations of y~'s are localized. Actually, at the leading 
order, the term proportional to%~-yf ) 0(-y 2 )-6(y 2 -y^) 6(y- -y 2 )-6(y^ -y 2 ) 6>(-yf) 
is the eikonal contribution in Feynman gauge to make the normal twist-2 quark distribution 
<p q /A(x) in Eq. flllD gauge invariant. Eq. ( p5|) is a good example to demonstrate that the 
double scattering does not give large nuclear size effect to the total inclusive cross section. 

On the other hand, from Eq. fl24|) , the double scattering contribution to the averaged 
transverse momentum square can gain large size effect 



A(q%) ~ J dq\q\ 



dQ 2 dq 2 J 
oc J dq^qq- $(<1t ~ ^t) ~ ^(Qt) 
~ k% . 



(26) 



Actually, k\ in Eq. fl2"BD needs to be integrated first. But Eq. ( |2"6| ) already demonstrates 
that A(g|>) is proportional k 2 -., which is the kick of the transverse momentum the parton 
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received from the additional scattering. The bigger the nuclear size, the bigger effective k\. 
As shown below, A(g|.) is proportional to the nuclear size. 

After demonstrating the physical picture of transverse momentum broadening in the 
Drell-Yan process, we are now ready to carry out the algebra. Working out the spinor trace 
for all three diagrams in Fig. |TJ, and combining three diagrams and drop the term localized 
in space, we obtain 



Mir) = J dq 2 T q 2 T (T X] / dx' (j)q/ h (x') J dx5(Q 2 - xx's) 

^^%--yr)0(-y 2 -)e<- + »- 

x d 2 y T r d2 ^ ikT . yT gi(fe 2 /( X >s))p+( v --v-) 



(27TY 

1 

2 



x - (p A \A + {y 2 ,0 T )^ q {0)j + tfj q {y~)A + {y^,y 1T )\p A ) 



x • [5{ql - 4) - *(<&)] (27) 

« a E / ^'^ A (x0 J dxT^x) 5{Q 2 - xx's) , (28) 

where o"o is defined in Eq. (|12|) and the four-parton correlation function is given by 



x ^ (PA|F a + (y 2 -)^(0)7 + ^(^)^ +a (yDbA) • (29) 

In deriving Eq. (|28|), we expend <5(gf. — A;^) at /c^ = 0, known as the collinear expansion, 
and keep only the first non- vanishing term which corresponds to the second order derivative 
term, 

5(q 2 T - k 2 T ) - 6(q 2 ) « S'(q 2 T )(-g^) . (30) 

We use the factor k^k^ in Eq. (^TJ) to convert the k^A + k^A + into field strength _p Q +_p/ 3 + 
by partial integration. Here, we work in Feynman gauge. The terms associated with other 
components of A p are suppressed by l/p + compared to those with A + , because of the 
requirement of the Lorentz boost invariance for the matrix elements ||. Th A {x) defined in 
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Eq. ( p9[) is the four-parton correlation function in a nucleus. The superscript (I) represents 
the initial state interaction, in order to distinguish from the similar four-parton correlation 
function defined in Eq. (p~a|) . More discussion on the relation between T q /A(x) and T^j A (x) 
will be given in Sec. IV. 

Combining Eqs. (|5[), flTT|), and (PSD, we obtain the nuclear enhancement of Drell-Yan 
(gf,) at the leading order in a s 

2 _ f^a s \ E q e\ldx'^ lh {x')T%{jlx')lx' 
{(lT) { 3 ) ' E q elJdx'<j ) g /h (x') ( j )q/A (T/x)/x' ■ 1 j 

III. JET BROADENING IN DEEPLY INELASTIC SCATTERING 

Consider the jet production in the deeply inelastic lepton-nucleus scattering, e{ki) + 
A(p) — > e(k2) +jet(l) +X. k\ and k 2 are the four momentum of the incoming and outgoing 
lepton respectively, and p is the momentum per nucleon for the nucleus with the atomic 
number A. With I being the observed jet momentum, we define the averaged jet transverse 
momentum square as 



*'^r^» AS- (32) 



? 2w i _ / ,,;> / > da e A I da e A 
dxsdQ^dtj. I dxsdQ 2 

where Xb = Q 2 / (2p ■ q), q = ki — k 2 is the momentum of the virtual photon, and Q 2 = —q 2 . 
The jet transverse momentum Ij- depends on our choice of the frame. We choose the Breit 
frame in the following calculation. Similar to Drell-Yan da/dQ 2 dq^, the jet transverse 
momentum spectrum, da/dxBdQ 2 dl^, is sensitive to the A 1 / 3 type nuclear size effect due 
to multiple scattering. On the other hand, the inclusive DIS cross section dojdxBdQ 2 = 
J dl\ do I dx BdQ 2 dl\ does not have the A 1 / 3 power enhancement. Instead, it has the much 
weaker A-dependence, such as the EMC effect and nuclear shadowing. To separate the 
multiple scattering contribution from the single scattering, we define the jet broadening as 

A(i 2 T ) = (i 2 T y A - (i 2 T y N . (33) 

Keeping only contribution from the double scattering, similar to Eq. @, we have 
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A(i 2 T ) -- /'-//:■/:■ dx d B a d Q 2dl 2 T I ■ (34) 

In the rest of this section, we derive the leading contribution to A(/^) by using the same 
technique used to derive the Drell-Yan A(g^) in last section. 
The general expression for the cross section in DIS is 



1 e 2 „ d 3 h 



= 7T 7V V W " U T^TB , (35) 



2s Q 4 (2tt) 3 2E. 
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with s = (p + ki) 2 . In Eq. (^3j), the leptonic tensor L^ u is given by the diagram in Fig. 0a, 

L Mi , = ^Tr(7 ■ ^7^7 • fc 2 7^) , (36) 

and W^ v is the hadronic tensor given by the diagram shown in Fig. |2|b. The leading order 
double scattering diagrams contributing to jet broadening are given in Fig. [3|. 

Feynman diagrams in Fig. ||| give the double scattering contribution to the hadronic 
tensor W^ u as, 

W£ = £ [dxdndx, I iklf^^MH^u^^l) , (37) 

with the matrix element 

rp(F) ( , \ f dy dyi dy2 d 2 y T 

1 v [X. X\. X 2 , for) — I ~, rrr 

V ' ' ' ; J 2tt 2tt 2tt (2tt) 2 

xe ixp+y- e «x 2 p + s/J e ik T -y T 

x ^ (P^I^(O) 7 + A+Gfe , Or) A+iyT, y 1T ) ^{y~) \p A ) , (38) 

where superscript (F) indicates the final-state double scattering. The matrix element 
is equal to the in Eq. ( |T6| ) if we commute the gluon fields with the quark fields P,|13 
In Eq. (|37|), H^ lu (x,Xi,X2,kT,p,q,l) is the corresponding partonic part. 
For the diagram shown in Figfja, the partonic part is given by 

jja — (J . JJ . \ . 

M " ^ V (xp + q) 2 + ie (xp + x 2 p + q) 2 — it 



x (2tt) 5{(xp + xxp + k T + q) 2 ) 5{l 2 T - k 2 T ) dl 2 T 
2tt 

{2p~ q y 



C ■ Hfiv • ~J^2. — ' 'H^r ^t) dl>T 



c-/ kr — 2q-kr, 1 1 

x^ar + ax-gB- T Q — ■ -, 39 

2p ■ q x — x B + ze x + x 2 —Xb — k 
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where C is the color factor for the diagrams in Fig. [| and given by C = 1/6. In Eq. (j39"D, 
the ^-function is from the phase space, and H^ v represents the spinor trace and is given by 

H^ v = (An) 2 a s a em e 2 Tr [7 • jry^ ■ (xp + x 2 p + q)Yl • h p l • {%P + q)lv\ P P Pa , (40) 

where p p p a is a result of our definition for the hadronic matrix element in Eq. (|38|) . 

Following the derivation of the Drell-Yan (q^), first, we carry out the integrations of the 
parton momentum fractions by using the 5-function and two poles in Eq. (j39l), 

H a = J dxdx 1 dx 2 e ixp+y ~ e ixip+ ^~ y ^ e iX2P+y * 5{l 2 T - k 2 T ) 

k\-2q-k T 1 1 



2p ■ q x — Xb + ie x + x 2 — XB — it 



x8(x + xi — xb — 

z,p ■ q x — Xb i" it 

= (2tt) 2 %f - y~) 6(y 2 ) e^+sr 5(l 2 T - k 2 T ) e\^) . (41) 

Similarly, we have the corresponding integrations for the interference diagram shown in 
Fig. 0b 

H b = Jdx dxi dx 2 e ixp+v ~ e iXlP+ ^- y ^ e iX2P+y * 5{l 2 T ) 

1 1 

xdix +xi — xb) ■ • , 2 n 1 

T — T d -4- If 1 ki,-2q-k T . 

x x B -t- tt x + x i — Xb ^ Vie 



i2n) 2 e(y 2 -y^9(y^-y-)e^ y -6(l 2 T )e\ 2p ' q J ' ' , (42) 



2p-q 

2 

krp—2q-krp 



p + (vi -y 2 ) 



and for the diagram in Fig. |3]c 



H c = J dxdx 1 dx 2 e ixp+y e ixip+ ^- y * } e iX2P+y ^ 5{l 2 T ) 
xd{x — xb) 



x + Xl -x B - T ~ 2p q ' q -ie x + x 2 -x B -te 



-{2nY9{y 2 )e{y^-y 2 )e^ p ^ S(l 2 T ) e V J . (43) 



Combining Eqs. (|4T|), (^) and ([£3|), we again have the same structure as that in Eq. 
Therefore, like the Drell-Yan case, we conclude that the double scattering does not result into 
any large nuclear size dependence to the inclusive DIS cross section. For the jet broadening 
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defined in Eq. (|34[), we drop the term proportional to — y~) 0(1/2) ~ 8(1/2 ~ Vi) 9(Vi ~ 
y~~) — 6(y± — y^) #(2/2"), which is localized as the single scattering. After carrying out the 
algebra similar to those following Eq. (p6|), we derive the jet broadening in DIS as 

where J2q sums over all quark and antiquark flavors. In Eq. (fill) , <f>q/A(%) is the normal twist- 
2 quark distribution inside a nucleus, and the four-parton correlation function, T q /a{xb) is 
defined in Eq. (|Ta]). 

From Eq. flH|), we conclude that based on the leading order calculation, we can uniquely 
predict the jet broadening in DIS, without any free parameter, if the four-parton correlation 
function, ^(ig), is measured in another experiment. Therefore, jet broadening in DIS 
can be used to test the QCD treatment of multiple scattering in nuclear environment. In 
addition, the jet broadening given in Eq. (H3) is directly proportional to ^(ig), in- 



dependence of A(iy) can provide immediate information on the functional form of the four- 
parton correlation function. 

IV. DISCUSSION AND CONCLUSIONS 

From Eqs. (|31~1) and flH]) , both jet broadening in DIS and the nuclear enhancement of the 
Drell-Yan (q^) are directly proportional to the four-parton (twist-4) correlation functions, 
which are defined in Eqs. flTa]) and (^), respectively. Since field operators are commuted 
on the light-cone PJI3H, the four-parton correlation functions, T q /^(x) and T^ A (x), are the 
same if the phase space integral are symmetric. With the model given in Eq. (Q), we predict 
the jet broadening in DIS as 

A</ 2 > = (^^j A 2 AW ; (45) 
and predict the nuclear enhancement of Drell-Yan transverse momentum square as 

A<<#) = I ^l) a 2 A 1/3 . (46) 
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which is the same as that obtained in Ref. [|J. It is a direct result of the model proposed 
by LQS |3|] that both jet broadening in DIS and nuclear enhancement in Drell-Yan A(g|,) 
have the very simple expressions, shown in Eqs. (^) and fl46|). In addition, Eqs. fl4"5|) and 
([46|) tell us that at the leading order, jet broadening in DIS and nuclear enhancement in 
Drell-Yan A(gf,) are the same, if the averaged initial-state gluon interactions is equal to the 
corresponding final-state gluon interactions, i.e., T q /^(x) = T^ A (x). 

From the simple expression in Eq. (|45|) , we conclude that at the leading order, the jet 
broadening in DIS has a strong scaling property, and it does not depend on beam energy, Q 2 
and xb- However, the x^-dependence needs to be modified if xb is smaller than 0.1. When 
x B < 0.1, the localized term, %f - y~) 6{y 2 ) - 6{y 2 - y±) %f - y-) - - y 2 ) 9{y 2 ), 
is no longer localized |ll],|14[]. Therefore, we have to keep this term for jet broadening 



calculation if xb is small. In addition, Q 2 -dependence may be modified because the four- 
parton correlation function T q /A{x) and the normal quark distribution <p q /A{x) can have 
different scaling violation. Of course, all dependence or whole conclusion could be modified 
due to possible different high order corrections. Nevertheless, we believe that experimental 
measurements of the jet broadening in DIS can provide valuable information on strength of 
multi-parton correlations and dynamics of multiple scattering. 

Similarly, from Eq. (f46[), we can also conclude that nuclear enhancement in the Drell-Yan 
(q 2 ^) has small dependence on beam energy and Q 2 of the lepton pair. However, data from 
Fermilab E772 and CERN NA10 |7],[| demonstrate some energy dependence. It signals that 
high order corrections for Drell-Yan (q 2 ^) is important P,10||, or the simple model of LQS 



for four-parton correlation functions is too simple. On the other hand, the weak energy 
dependence from the data shows that the leading order calculation given here can be useful, 
and QCD treatment of multiple scattering can be eventually tested and understood. 

Use Eq. ( f46|) and data from E772 and NA10 on nuclear enhancement of Drell-Yan (q^), 
we estimate the value of A 2 as 

X 2 DY w O.OlGeV 2 , (47) 
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which is at least a factor of five smaller than \%_ }et ~ 0.05 ~ 0.1 GeV 2 , estimated by LQS 
from momentum imbalance of di-jet data ||. 

Since jet broadening in DIS and the momentum imbalance in di-jet are both due to the 
final-state multiple scattering, we will use the value of A? _ jct to predict the jet broadening as 

A(^) « (0.66 ~ 1.31) a s A 1 ^ . (48) 

On the other hand, with a smaller value of A^ Y in Eq. (£|7]), we will predict a much smaller 
jet broadening ~ 0.13 a s A 1 / 3 . Direct experimental measurement on jet broadening can 
certainly provide an independent test of our QCD treatment of multiple scattering. Early 
data on jet production from Fermilab E665 can be used to calculate the jet broadening. 



Future experiments at HERA with a heavy ion beam |TTJ should be able to provide much 
more information on dynamics of parton correlations. 

In summary, the jet broadening in DIS provides an independent measurement of four- 
parton correlation functions, and tests of QCD dynamics beyond simple parton model. At 
the leading order, jet broadening is directly proportional to the four-parton correlation 
functions. The size of jet broadening provides a direct measurement of the size of the 
four-parton correlation functions. In addition, by varying the xb, we can gain valuable 
information on x-dependence of the four-parton correlation functions. Information on in- 
dependence of the jet broadening, A(/y) in Eq. (|44]) , can provide a direct measurement of 
the functional form of T q / A (xB, A) defined in Eq. (|Ta|). 
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FIGURES 

FIG. 1. Lowest order double scattering contribution to the nuclear enhancement of Drell-Yan 
((/r); (a)symmetric diagram; (b) and (c): interference diagrams. 

FIG. 2. Diagrams for DIS: (a) Diagram representing L^; (b) Diagram representing W^ u . 

FIG. 3. Lowest order double scattering contribution to jet broadening: (a) symmetric diagram; 
(b) and (c): interference diagrams. 
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